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Abstract 

We summarize some of the recent developments which link certain prob- 
lems in combinatorial theory related to random growth to random matrix 
theory. 

2000 Mathematics Subject Classification: 60C05. 

Keywords and Phrases: Random matrices, Toeplitz determinants, Deter- 
minantal processes. Random Growth, Random permutations. 



1. Introduction 

Let CT be a permutation from Sn- We say that cr(ii), . . . , aiim), ii < ■ ■ ■ < im, 
is an increasing subsequence of a if cr(ji) < • • • < a{im)- The number m is the 
length of the subsequence. The length of the longest increasing subsequence in a 
is denoted by In (a). If we pick a from Sn uniformly at random £Ar(cr) becomes a 
random variable. Ulam's problem, [221 j is the study of the asymptotic properties 
as iV ^ oo of this random variable in particular its mean. It turns out that 
there is a surprisingly rich mathematical structure around this problem as we hope 
will be clear from the presentation below. It has been known for some time that 
E[£jv] ~ 2\/]V as iV ^ oo, HD, [E]. We refer to for some background to the 
problem. A Poissonized version of the problem can be obtained by letting N be an 
independent Poisson random variable with mean a. This gives a random variable 
L{a) with distribution 

OP —a N 

P[L{a) <n]^Yl ^^^[^^ ^ (1-1) 

N=0 

Since V[£n < n] is a decreasing function of A^, asymptotics of the left hand side 
of (1.1) can be used to obtain asymptotics oiP[£N < n] (de-Poissonization). 
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The random variable L{a) can be reahzed geometrically using Hammersley's 
picture, j8|. Consider a Poisson process with intensity 1 in the square [0,7]^, 7 = 
y/a. An up/right path is a sequence of Poisson points (xi, yi), . . . , (a^m, J/m) in 
the square such that Xi < Xi+i and yi < Ui+i, i — 1,...,to — 1. The maximal 
number of points in an up/right path has the same distribution as L(a). A sequence 
of points realizing this maximum is called a maximal path. It is expected from 
heuristic arguments, see below, that the standard deviation of L{a) should be of 
order 7^/'^ = a^^^. The proof that this is true, [H], and that we can also understand 
the law of the fluctuations is the main recent result that will be discussed below. 
Also, the deviations of a maximal path from the diagonal x — y should be of order 
7^/'^. This last statement is proved in 

A generalization of the random variable L{q) can be defined in the following 
way. Let wii,]), S Z^, be independent geometric random variables with 

parameter q. An up/right path tt from (1,1) to {M,N) is a sequence (1,1) = 
(n, ii), («2,j2), («m,im) = (M, iV) , TO = M+iV- 1 , such that either v+i-V = 1 
and jr+i = jr, or v+i = v and j^+i — jr = 1- Set 

G(M,iV)= max V w{tj), (1.2) 

7r ^ — ^ 

where the maximum is taken over all up/right paths tt from (1, 1) to {M,N). Al- 
ternatively, we can define G{M, N) recursively by 

G{M, N) = max(G'(Af - 1, N),G{M, - 1)) + w{M, N). (1.3) 

Some thought shows that if we let q = a/N"^ then G{N,N) converges in distri- 
bution to L{a) e.s N ^ oo, |12j, so we can view G{N,N) as a generalization of 
L{a). We can think of (|1.2|l as a directed last-passage site percolation problem. 
Since all paths tt have the same length, if w{i,j) were a bounded random variable 
we could relate ()1.2|l to the corresponding first-passage site percolation problem, 
with a min instead of a max. in (|1.2|l . The random variable G{M,N) connects 
with many different problems, a corner growth model, zero-temperature directed 
polymers, totally asymmetric simple exclusion processes and domino tilings of the 
Aztec diamond, see ^01: and references therein. It is also related to another 
growth model, the (discrete) polynuclear growth (PNG) model, |15| . |2H defined as 
follows. Let h{x,t) S N denote the height above a; G Z at time t gN. The growth 
model is defined by the recursion 

h{x, t + 1) = max(/i(x - 1, t),h{x, t), h{x + 1, t)) + a{x, t), (1.4) 

where a{x,t), {x,t) S Z x N, are independent random variables. If we assume 
that a{x,t) = whenever x — t is even, and that the distribution of a{x,t) is 
geometric with parameter q, then setting w{i,j) — a{i — j,i + j — 1), we obtain 
G{i,j) = h{i — j,i + j — 1). The growth model (|1.4|) has some relation to the so 
called Kardar-Parisi-Zhang (KPZ) equation, 15 , and is expected to fall within the 
so called KPZ-universality class. The exponents 1/3 and 2/3 discussed above are 
the conjectured exponents for 1 + 1-dimcnsional growth models in this class. 
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2. Orthogonal polynomial ensembles 

Consider a probability (density) on 17^, Q = M, Z, N or {0, 1, . . . , M} of the 



form 

1 

un{x) = —ANixf TT w{xj), (2.5) 

J = l 

where An(x) = Y[i<i<j<Ni^i Vandermonde determinant, w{x) is some 

non- negative weight function on Q and is a normaUzation constant. We call 
such a probability an orthogonal polynomial ensemble. We can think of this as a 
finite point process on f2. Let dfj, be Lebesgue or counting measure on f2 and let 
Pn{x) = Kri^;" + . . . be the normalized orthogonal polynomials with respect to the 
measure w{x)d^{x) on fi. The correlation functions Pm.N{xi, ■ ■ ■ ,Xm) of the point 
process are given by determinants, we have a so called determinantal point process, 
m. In fact, 

iV! f 

Pm.N{xi, ■ ■ ■ ,Xm) = — r- / UN{x)dfJ.{x.m+l) ■ ■ -dnixN) (2.6) 

= det{KN{xi,Xj))i<ij<„i, 
where the kernel Kn is given by 

KN-lPN{x)pN-l{y) - PN-l{x)pN{y) , , . I ^.l/2 /o7\ 

KN{x,y)^ {w{x)w{y)) ' . (2.7) 

KN X - y 

A computation shows that for bounded f : —> C, 

N N , ^ k 



ntli^ + f{x,))]^J2h /^n/(^.)det(X^(a;.,x,))i<.,<ferfV(a:) (2.8) 
j=i k=o ■ j=l 

= det(/ + fKN)mn,dt,), 



where the last determinant is the Fredholm determinant of the integral operator 
on L^{fl, dfi) with kernel f{x)KN{x, y). In particular, we can compute hole or gap 
probabilities, e.g. the probability of having no particle in an interval / C 57 by taking 
/ = —xii minus the characteristic function of the interval /. If Xmax — maxccj 
denotes the position of the rightmost particle it follows that 

IP[a;max < a] = det(/ - i^Ar)L2((„ ^o) j^^). (2.9) 

As N oo we can obtain limiting determinantal processes on R or Z with 
kernel i.e. the probability (density) of finding particles at xi, . . . ,Xm is given 
by det{K{xi,Xm))i<i,j<m- We will be interested in the limit process around the 
rightmost particle. This is typically given by the Airy kernel, a;, y € M, 

^^^^y^^mx)M'iy)-M'ix)Miy)_ 
x-y 
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This hmit, the Airy point process has a rightmost particle almost surely and its 
position has the distribution function 



known as the Tracy- Widom distribution, |26j . Hole probabilities as functions of 
the endpoints of the intervals satisfy systems of differential equations, |22], For 
example, we have 



where u solves the Painleve-II equation u" ~ xu + 2u'^ with boundary condition 
u(x) ^ Ai{x) as X — > oo. 

An example of a measure of the form (|2.5|l comes from the Gaussian Unitary 
Ensemble [GUE) of random matrices. The GUE is a Gaussian measure on the 
space Ti-N = of all x iV Hermitian matrices. It is defined by (i/XGUE.Ar(M) = 
ZJ^^ exp{—tr M^)dM , where dAI is the Lebesgue measure on TCm and is a 
normalization constant. The corresponding eigenvalue measure has the form (|2.5|l 
with w{x) = exp(— a;^) and = M, 17 . Hence the p„:s are multiples of the 
ordinary Hcrmite polynomials. The largest eigenvalue Xmax will lie around \/2N . 
This is related to the fact that the largest zero of pn lies around V2iV. The local 
asymptotics oipN{x) exp(— x^/2) around this point, x — ^/2N + ^/7V^/^-\/2, is given 
by the Airy function, Ai (^). This asymptotics, some estimates, (|2.9ll and H2.11|l give 
the following result, 



as A^ ^ oo. 

3. Some theorems 

The previous section may seem unrelated to the first but as the next theorems 
will show the problem of understanding the distribution of L{a) and G{M, N) fits 
nicely into the machinery of sect. 2. 

Theorem 3.1. ilO|. Take n ^N, M > N and w{x) = (^'^"f in 
Then G{M,N) is distributed exactly as Xmax- 

The corresponding orthogonal polynomials are the Meixner polynomials, a 
classical family of discrete orthogonal polynomials, and we refer to the measure 
obtained as the Meixner ensemble. It is an example of a discrete orthogonal polyno- 
mial ensemble, |12j . By computing the appropriate Airy asymptotics of the Meixner 
polynomials we can use ()2.9|l to prove the next theorem. 



Theorem 3.2. |Tn|. Let-f > I be fixed and set uj{-f ,q) = {l-q)-^{l+^Y -I 
and a(7, q) ^ {I - qY^qhY'^^ + ^f'^{l + ^f'^- Then, 



F2(0 = det(/-^)i2(5,^), 



(2.11) 




(2.12) 




(2.13) 



(3.14) 
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Thus G'([7A''], A'') fluctuates Hke the largest eigenvalue of a GUE matrix. 

As discussed above by setting q = a/N'^ we can obtain L{a) as a limit of 
G{N,N) a.s N ^ oo. By taking this limit in theorem 3.1, and using the fact that 
the measure has dcterminantal correlation functions, we see that L{a) behaves like 
the rightmost particle in a determinantal point process on Z given by the discrete 
Bessel kernel, 0], 

B {x,y)^^/a , (3.15 

x-y 

x,y G Z. This gives 

P[L(a) <n] =det(/-B")p({„,„+i^...}). (3.16) 

Once we have this formula we see that all we need is the classical asymptotic formula 
Q"'^/^ J2^_i_^Q,i/6 (2y^ —>■ Ai (^) as a — > cxD uniformly in compact intervals, and some 
estimates of the Bessel functions in order to get a limit theorem for L{a): 
Theorem 3.3. 0. As a oo, 

n ^ ^1/6 <i\^ my (3.17) 

Note the similarity with (|2.13|) . just replace N by ^fa. This result was first 
proved in [3] by another method, see below. De-poissonizing we get a limit theorem 
for i?Ar(7r), see 



4. Rewriting Toeplitz determinants 

The Toeplitz determinant of order n with generating function / G L^(T) is 
defined by 

A. (/) =det (/,_,) (4.18) 

where — (27r)~^ j(e^^)e^^^^ dQ are the complex Fourier coefficients of /. Con- 
sider the generating function 

M 



f{z)^ll{l + ^)il + hz), (4.19) 
e=i ^ 

where ai,b£ are complex numbers. The elementary symmetric polynomial em(a), 
a = (ai, . . . , gm) is defined by njii(l + o.jz) = J2\,n\<oo em(a)z'"- A straightfor- 
ward computation shows that when / is given by H4.19|) then 

oo 

fi-j = ^ e„i-j{a)era-k{b). 



m=0 
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Insert this into the definition H4.18|l and use the Heine identity, 

^ / det{<j)i{xj))l<^J<ndet{'^pi{xj))l<i^J<nd"^i{x) (4.20) 

= det( / (l}i{x)i)j{x)d^i{x))i<ij<ri, 
Jn 

to see that 

Dn{f)= ^ det(em,_j(a))i<ij<„ det(e„i,_j(&))i<ij<„. (4.21) 

mi >m2 >'">m^>0 

Here we have removed the n\ by ordering the variables. These determinants are 
again symmetric polynomials, the so called Schur polynomials. Let A = (Ai, A2, . . . ) 
be a partition and let A' = (A^ , A2, . . . ) be the conjugate partition, |23. Set m.^ = 
X[-\-n — i,i = \,...,n and AJ = if i > n, so that A' has at most n parts, ^(A') < n, 
which means that Ai < n. Then the Schur polynomial s\{a) is given by 

sx{a) = det(eA^_j+j(a))i<,jj<„ = det(emi-i(a))i<ij<„, (4.22) 
the Jacobi-Trudi identity. Hence, 

A.(/)= •^A(a)sA(fe), (4.23) 

A ; Ai <n 

and we have derived Gessel's formula, "T*. If we let n 00 in the right hand side 
we obtain Jli ~ ciibj)^^ by the Cauchy identity, In the case when all 

Ui, bj G [0, 1], sx{a)sx{b) > 0, and we can think of 

M 

Y[{l-a,bj)sx{a)sx{b) (4.24) 

as a probability measure on all partitions A with at most n parts, the Schur mea- 
sure, |19j . In this formula we can insert the combinatorial definition of the Schur 
polynomial, j23| . 

s.ia)= E ar^"'..-r^"\ (4.25) 

T:sh(T)=A 

where the sum is over all semi-standard Young tableaux T, |23j . with shape A, and 
mi{T) is the number of i:s in T. 

A connection with the random variables in section 1 is now provided by the 
Robinson-Schensted-Knuth (RSK) correspondence, This correspondence maps 
an M X M integer matrix to a pair of semi-standard Young tableaux (T, S) with 
entries from {1, 2, . . . , M}. If we let the random variables w{i,j) be independent 
geometric with parameter aibj then the RSK-correspondence maps the measure we 
get on the integer matrix {w{i, j))i<ij<M to the Schur measure 14.24|l . Also, the 
RSK-correspondence is such that G{M, M) — Ai, the length of the first row. If we 
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put aj ~ for N < j < M and a.i = hi = y/g for 1 < i < TV in the Schur measure 
and set Xj = Xj+N — j, 1 < j < N, we obtain the result in Theorem 3.1. 

In the Hmit M ^ N ^ oo, g = a/N'^, in which case G{N,N) converges to 
L(a), the Schur measure converges to the so called Plancherel measure on partitions, 
|30|. |12|. In the variables Xi — i this measure is a dcterminantal point process on 
Z given by the kernel B", H3.15|l . This result was obtained independently in 0], 
which also gives a descrption in terms of different coordinates. See also JHl for a 
more direct geometric relation between GUE and the Plancherel measure. In this 
limit the Toeplitz determinant formula (|4.23|) gives 



This variant of Gessel's formula was the starting point for the original proof of 
Theorem 3.3 in The right hand side of 14.26|l can be expressed in terms of the 
leading coefficients of the orthogonal polynomials on T with respect to the weight 
exp(2Y^cos 9). These orthogonal polynomials in turn can be obtained as a solution 
to a matrix- valued Riemann-Hilbert problem (RHP), and the asymptotics of this 
RHP as a — > cx) can be analyzes using the powerful asymptotic techniques developed 
by Deift and Zhou, This approach leads to the formula l|2.12|l for the limiting 
distribution. 

Write / = exp((7) and insert the definition of the Fourier coefficients into 
the definition H4.18I) . By the Heine identity we obtain an integral formula for the 
Toeplitz determinant, 



In the last integral dU denotes normalized Haar measure on the unitary group U (n) 
and the identity is the Weyl integration formula. The limit of H4.26|l as a ^ oo is 
then a so called double scaling limit in a unitary matrix model, |20| . The formula 
(|4.26|l can also be obtained by considering the integral over the unitary group, see 



Another way to obtain the Schur measure is via families of non-intersecting 
paths which result from a multi-layer PNG model, The determinants in the 

measure then come from the Karlin-McGregor theorem or the Lindstrom-Gessel- 
Viennot method. 

5. A curiosity 

Non-intersecting paths can also be used to describe certain tilings, e.g. domino 
tilings and tilings of a hexagon by rhombi. By looking at intersections with appro- 
priate lines one can obtain discrete orthogonal polynomial ensembles. In the case 
of tilings of a hexagon by rhombi, which correspond to boxed planar partitions, [S], 



P[L(a) < n] = e-"D„(e2v^™^''). 



(4.26) 




(4.27) 



22 . 
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the Hahn ensemble, i.e. H2.5|l with = {0, . . . , M} and a weight giving the Hahn 
polynomials, is obtained, J^I- The computation leading to this result also gives 
a proof of the classical MacMahon formula, [221 j for the number of boxed planar 
partitions in an abc cube, i.e. the number of rhombus tilings of an a5c-hexagon. In 
terms of Schur polynomials the result is 

E v(nv(i-)-nnn ;;^:;^:; . (5.28) 

fj.;e{ii)<c i=lj=lk=l ■' 

where the right hand side is MacMahon's formula. (Here 1" means (!,...,!) with 
a components.) Comparing this formula with the formula 14.23|) we find 

i-ir'D^ifia - e--) rid - e-)) = n ri n (5-29) 

e=i 1=1 i=ij=ik=i ■' 

It has been conjectured by Keating and Snaith, that the following result should 
hold for the moments of Riemann's (^-function on the critical line, 

1 1 



TT^^?^ / IC(l/2 + ^i)l''rfi = /cuE(fc)a(fc), (5.30) 
T^oo (logT)'' T Jq 

where a{k) is a constant depending on the primes, 

-IP ^ ' I 

/cuE(fc) = lim — / \ZiU,e)f''dU = n r^TTT (^.31) 

and Z{U, 9) — det(/ — C/e^*^) is the characteristic polynomial of the unitary matrix 
U. If we take a = 6 = fc in H5.29|l and use (|4.27(l we find 

as computed in ' 14| by different methods. Letting n — > cxd we obtain the last expres- 
sion in (|5.31|l . Hence, we see that the formula (|5.32|l has a curious combinatorial 
interpretation via MacMahon's formula. 
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